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ABSTRACT 

A star that wanders too close to a massive black hole (BH) is shredded by the BH's 
tidal gravity. Stellar gas falls back to the BH at a rate initially exceeding the Eddington 
rate, releasing a flare of energy. In anticipation of upcoming transient surveys, we 
predict the light curves and spectra of tidal flares as a function of time, highlighting the 
unique signatures of tidal flares at optical and near- infrared wavelengths. A reasonable 
fraction of the gas initially bound to the BH is likely blown away when the fallback 
rate is super-Eddington at early times. This outflow produces an optical luminosity 
comparable to that of a supernova; such events have durations of ~ 10 days and may 
have been missed in supernova searches that exclude the nuclear regions of galaxies. 
When the fallback rate subsides below Eddington, the gas accretes onto the BH via 
a thin disk whose emission peaks in the UV to soft X-rays. Some of this emission 
is reprocessed by the unbound stellar debris, producing a spectrum of very broad 
emission lines (with no corresponding narrow forbidden lines). These lines are the 
strongest for BHs with Mbh ~ 10 5 — 10 6 Mq and thus optical surveys are particularly 
sensitive to the lowest mass BHs in galactic nuclei. Calibrating our models to ROSAT 
and GALEX observations, we predict detection rates for Pan-STARRS, PTF, and 
LSST and highlight some of the observational challenges associated with studying 
tidal disruption events in the optical. Upcoming surveys such as Pan-STARRS should 
detect at least several events per year, and may detect many more if current models of 
outflows during super-Eddington accretion are reasonably accurate. These surveys will 
significantly improve our knowledge of stellar dynamics in galactic nuclei, the physics 
of super-Eddington accretion, the demography of intermediate mass BHs, and the role 
of tidal disruption in the growth of massive BHs. 
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1 INTRODUCTION 

Stellar orbits in the center of a galaxy are not static, and 
sometimes stars walk into trouble. If an unlucky star passes 
within Rt ~ R+^Mbh/M*) 1 / 3 of the galaxy's central black 
hole (BH), the BH's tidal gravity exceeds the star's self- 
gravity, and the star is disrupted. For BHs with Mbh ^ 
10 s Mq, the disruption of a solar- type star occurs outside 
the horizon and is likely acco mpanied by a week- to year- 
long electromagnetic flare fe.g.. lReeslll988l h 

Gravitational interactions between stars ensure that 
all supermassive BHs tidal ly disrupt nearby stars (e.g., 
iMagorrian fc Tremainej[l999l ). The scattering process might 
be accelerated by resonan t relaxation very close to the BH 
(|Rauch fc Tremainel 1 19961 ). or interactions with "massive 
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perturbers" like a massive accretion disk dZhao et alj|2002h 
or giant molecular clouds l|Perets et al.ll2007f ). In addition, 
the galactic potential may be triaxial so stars need not be 
scattered at all: they may simply follow their chaotic orbits 
down to ~ Rt (|Merritt fc Poonl 120041 ') . Given these uncer- 
tainties, predictions for the timescale between tidal disrup- 
tions in a given galaxy range from 10 3 to 10 6 years. The rate 
remains uncertain, but tidal disruption must occur. 

Indeed, a handful of candidate events have been de- 
tected to date. The accreting stellar debris is expected to 
emit blackbody radiation from very close to the BH, so 
X-ray and UV observations probe the bulk of the emis- 
sion. Several candidate tidal disruption eve nts were discov - 
ered in the ROSAT All-Sky Survey (see iKomossal l2002h 
and the XMM-Newton Slew Survey (|Esquei et alj 120071 ): 
the GALEX Deep Imaging Survey has so fa r yielded three 
candidates (jGezari et al.l |2006| . 120081 . [2009). For ROSAT, 
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these d etections are consist ent with a rate ~ 10 -5 yr~ per 
galaxy (|Donlev et al.l [20021) . but the data are sparse. How- 
ever, we are entering a new era of transient surveys: in the 
optical, surve ys like Pa n-STARRS (PS1, then all four tele- 
scop esj_(e ; g_ li jMagnie3 [2007]), the Palomar Transient Fac- 
tory (|Rau et alj|2009l ). and later the Large Synoptic Survey 
Telescope will have fast cadence, wide fields of view, and 
unprecedented sensitivity. Wide-field transient surveys with 
rapid cadence are also planned at other wavelengths, includ- 
ing the radio (e.g., LOFAR and the ATA), near-infrared 
(e.g., SASIR), and hard X-rays (e.g., EXIST). How many 
tidal flares these surveys find depends on their luminosity 
and spectra as a function of time. 

In this paper, we predict the light curves and spectra 
of tidal disruption events as a fun ction of tim e. Since the 
early work on tidal disruption (e.g.. lReeslll988t ). it has been 
well-appreciated that the bulk of the emission occurs in the 
UV and soft X-rays, with a possible extension to harder X- 
raysQ Taking into account only this emission, optical wave- 
lengths are not the most promising for detecting tidal flares, 
because the blackbody temperature of the inner accretion 
disk is ~ 3 x 10 K. We show, however, that there are two 
additional sources of optical emission that likely dominate 
the optical flux in many, though not all, cases: (1) emission 
produced by a super-Eddington outflow at early times and 
(2) emission produced by the irradiation and photoioniza- 
tion of unbound stellar debris (see also the earlier work of 
Bogd anovic et ahll2004l ). In the latter case, much of the opti- 
cal emission is in the form of very broad emission lines, while 
in the former it is primarily continuum (although some lines 
may also be present). Throughout this paper, we typically 
discuss these two sources of emission separately, largely be- 
cause the physics of the photoionized stellar debris is more 
secure than that of the super-Eddington outflows. 

The remainder of this paper is organized as follows. 
In ^2.11 12.21 and [3] we describe our models for the polar 
super-Eddington outflow, accretion disk and the equatorial 
unbound material, respectively; then in Sj4]we calculate the 
luminosity and spectral signatures of tidal disruption events. 
We predict detection rates in 33 and summarize and discuss 
our results in ^6] § i|5.2l and [6] include a discussion of our 
models in the context of ROSAT and GALEX observations 
of tidal flare candidates. 



2 THE INITIALLY BOUND MATERIAL 

We consider a star approaching the BH on a parabolic orbit 
with pericenter distance R p ^ Rt- Once the star reaches the 
vicinity of the tidal radius (-Rt), the tidal gravity stretches 
it radially and compresses it vertically and azimuthally. The 
acceleration is a ~ (GMbh/Rp){R*/R p ) and acts for a dy- 
namical time t p ~ (GMbh/Rp)" 1 ^ 2 near pericenter, result- 
ing in velocity perturbations Av p ~ at p ~ Vi,(RT / R P ) 3 , 

1 Such a hard X-ray component may be detectable with up- 
coming all-sky X-ray surveys like the proposed EXIST mission 
llGrindlay||2003) . However, we choose not to include predictions 
for hard X-rays in our calculations: one could draw analogy to the 
hard X-ray power-law tail observed in AGN spectra, but since the 
origin of this feature is uncertain, it is difficult to make firm the- 
oretical predictions for tidal disruption events. 



where R+ is the star's radius and -u* is the star's escape ve- 
locity. The change in velocity Av p is smaller than the star's 
orbital velocity at pericenter, v p = (2GMbh/ Rp) 1 ^ 2 , by a 
factor of Rir/Rp. 

Because Av p is at least as large as the sound speed 
inside the star, the stellar gas may shock vertically and az- 
imut hally (e.g., Br assart fc Luminetll200§ : iGuillochon et al.l 
2008). Once the shredded star passes through pericenter, 
the compression subsides and the star re-expands, cooling 
adiabatically; thermal pressure becomes negligible and the 
particles travel away from the BH ballistically. We assume 
that the particle trajectories become ballistic when the star 
passes through pericenter. At that time, the particles have 
perturbed azimuthal, vertical, and radial velocities ~ Av p . 

The particles have a ran ge in speci fi c ene r gy £ ~ 
±3(G , M BH /-Rp)(^/-Rp) (e.g.. lLacv et all Il982l : iLi et all 
2002), due to their relative locations in the BH's potential 
well and differences in their azimuthal speeds. Initially, ap- 
proximate ly half of the stellar mass is bound and half is 
unbound (|Lacv et al.lll982l ; lEvans fc Kochaneklll989h . After 
a time 

2n (R P \ 3/ \ m 

~ 20M 6 5/2 < 3Rsr -3/2 mm] 

the most bound material returns to pericenter. (Here we 
have defined M 6 = M BH /1O 6 M , Rp,m. s = R P /SR S , and 
r* = R+/Rq.) Less bound gas follows, at a rate 

-Mfallback ~ 77 ~ I 7 (2) 

o tfallback V tfallback / 

l|Reeslll988l : iPhinnevI fl989h . There will be some deviations 
from this canonical t~ 5 ^ 3 scaling at e arly times, depend- 
ing on the precise structure of the star (|Lodato et al.ll2009l : 
iRamirez-Ruiz fc Rosswo"gll2009l ) , but we use equation ([2J for 
simplicity. As matter returns to pericenter, it shocks on it- 
self, conver ting most of its bulk orbital energy to thermal 
energy (see lKochanek|[l994l ). The viscous time is typically 
shorter than the fallback time, so at least some of the matter 
begins to accrete. 

For Mbh <^ few x 10 7 Mq, the mass fallback rate pre- 
dicted by equations (fTJ & @ can be much greater than the 
Eddington rate MEdd for a period of weeks to years; here 
M E dd = 10I/Edd/c 2 , I/Edd is the Eddington luminosity, and 
0.1 is the assumed efficiency of converting accretion power to 
luminosity. The fallback rate only falls below the Eddington 
rate at a time 

j f. -| „2/5 D 6/5 3/5 -3/5 ro \ 

t E dd~0.1M 6 ' Jip' 3Hs m/ r* yr, (3) 

where m* = M+/Mq. While the fallback rate is super- 
Eddington, the stellar gas returning to pericenter is so dense 
that it cannot radiate and cool. In particular, the time for 
photons to diffuse out of the gas is longer than both the in- 
flow time in the disk and the dynamical time characteristic 
of an outflow. The gas is likely to form an adve ctive accretion 
disk accompanied by powerful outflows (e.g.. lohsuga et "ail 
12005! ) . although the relative importance of accretion and 
outflows in this phase is somewhat uncertain (see Later, 
when A/f a nb a ck < A^Edd (t > *Edd), the outflows subside, and 
the accretion disk can radiatively cool and becomes thin. 
In Q2.1\ we describe our model for the super-Eddington 
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outflows, and in £12.21 we describe our model for the accretion 
disk. We discuss uncertainties in these models in 36] 



2.1 Super-Eddington Outflows 

When the fallback rate to pericenter is super-Eddington, 
radiation produced by the shock and by viscous stresses 
in the rotating disk is trapped by electron scattering. By 
energy conservation, this material is initially all bound to 
the BH, but it is only weakly bound because the radia- 
tion cannot escape and because the material originated on 
highly eccentric orbits. Some fraction of the retu rning gas is 
thus likely unbound (see, e.g., the simulations of lAval et al.l 
2000), with energy being conserved as other gas accretes in- 
ward l|Blandford fc Begermanlll999l ). If the outflow's cover- 
ing fraction is high, most of the radiated power will be emit- 
ted fro m the outflow's photo sphere, which can be far outside 
~ R p jLoeb fc Ulmerlll997l ). We now estimate the prop er- 
ties of this outflowing gas (see iRossi fc Beg elman 20QSJ for 
related estimates in the context of short-duration gamma- 
ray bursts) . 

In our simplified scenario, stellar debris falls back at 
close to the escape velocity and shocks at the launch radius 
Rl ~ 2_R P , converting bulk kinetic energy to radiation: 



rr4 1 2 

" 2 /3fallback ' L ' Ucsc ' L ' 



(4) 



where Tl is the temperature at Rl, £>fallback,L ~ 
Mf a iiback/(47r_RL' !; csc,L) is the density of gas at _Rl, and 
^csc.L = (2GM bh /Rl) 1/2 . Outflowing gas is launched from 
Rl, at a rate 



M 



out 



/outMj 



fallback 



and with terminal velocity 

Vwind = fvVesc(Rh)- 



(•>) 



(6) 



We approximate the outflow's geometry as spherical, with a 
density profile 



p(r) 



47rr- 2 u wind 



(7) 



inside the outflow where r <J R c d gc = Vwtndi; the density falls 
quickly to zero at ~ -Rcdgc- We define the trapping radius 
-Rtrap via Rtiapp{Rtiap)n B ~ c/u w ind (where k, b is the opac- 
ity due to electron scattering): inside Rtrap, the gas is too 
optically thick for photons to escape and so the outflowing 
gas expands adiabatically. Because the outflow remains sup- 
ported by radiation pressure, T tx p 1 ^ 3 . The photosphere of 
the outflow R pri is where R! prl p(-Rph)K s ~ 1. Because u W md is 
likely not much smaller than c, -Rtrap ~ R P h; we thus neglect 
any deviations from adiabaticity between -Rtrap and R p h. 

At the earliest times for small A/bh and R p , the fallback 
rate can be so large and the density so high that the edge of 
the outflow limits the location of the photosphere to be just 
inside R cd gc. In that case, the density of the photosphere 
is still given by p(R p h) ~ (^s-Redgc) -1 ; lacking a detailed 
model, we assume that the photospheric gas near the edge 
is on the same adiabat as the rest of the gas, so that 



T, 



P h 



P(-Rph) 



Pfallback.L 



1/3 



3 x 10 4 /„" 1/3 M 6 1/36 i? 



■1/8 
p,3Rs 



-1/12 1/12 



day / 



-7/36 



K. 



Note that the photospheric temperature during the edge- 
dominated phase is essentially independent of all parameters 
of the disruption (e.g., Mbh, -Rp, etc.), and is only a weak 
function of time. The total luminosity during this phase 
grows as L tx t 11 ' 9 while the luminosity on the Rayleigh- 
Jeans tail increases even more rapidly, vL v tx t 65/36 . After 
a time 



^cdg 



1 f 3/8 ,-3/4 
^ J out Jv 



M, 



5/8 „9/8 



fl p',3R s m 



3/8-3/8 



day, 



(9) 



the density falls sufficiently that the photosphere lies well 
inside -R d gc ; the photosphere's radius is then 



-Rph ~ 4/out/u 1 



Mn 



M^dd 

and its temperature is 



(10) 



T, 



P h 



2xio 5 /o" u i / 3 /y 3 



M, 



M E 



-5/12 



The adiabatically expanding outflow preserves the photon 
distribution function generated in the shock and accretion 
disk close to the BH. Estimates indicate that this gas is likely 
to be close to thermal equilibrium and thus we assume that 
the escaping photons have a blackbody spectrum 



vL v ~ 4tt -Rph^-Bi/(Tph) • 



(12) 



When the photosphere lies inside the edge of the outflow 
(i.e., t > tcdgc so Rph < u w ind£), equations ((10} and (JTTJ) 
imply that the total luminosity of the outflow is 

L ~ 10 44 /o^ 3 / 1 T 2/3 M 6 11 / 9 <^ s my3 r -V3 (JL)~ B/8 ergs^.ClS) 

The total luminosity of the outflow is thus of order the Ed- 
dington luminosity: see Figure |4j discussed in 34.11 Note 
that the total luminosity decreases for lower outflow rates, 
L tx /^ t 3 , because the photosphere's surface area is smaller. 
The luminosity on the Rayleigh- Jeans tail (generally appro- 
priate for optical and near-infrared wavelengths) declines 
even faster for lower M ou t, scaling as 



vL v tx fUtfv 



,/3 



(14) 



(8) 



These relations only apply if _R p h ;> J?l because otherwise 
the outflow is optically thin; we impose this lower limit to 
-Rph in our numerical solutions described later. 



2.2 The Accretion Disk 

We now consider the bound stellar debris that accretes onto 
the BH. After shocking at pericenter, this gas circularizes 
and viscously drifts inward, forming an accretion disk. The 
disk extends from ~ 2R P down to the last stable orbit, 
Rlso- We expect the viscous time in the disk to be sub- 
stantially shorter than the fallback time for at least a few 
years Csee lUlmer![l999l for the case of R p = Rt, assuming a 
thick disk) , and check this expectation at the end of H2.2I we 
thus assume that accretion during this period proceeds at 
M ~ (1— /out)MfaUback- During the super-Eddington phase, 
the time for photons to diffuse out of the disk is longer than 
the viscous time, and so the disk is thick and advective. In 
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contrast, at later times when Mf a iiback <J Afedd, the disk 
is thin and can cool by radiative diffusion. We derive an 
a nalytic "slim disk" mode l (similar to the numerical work 
oflAb ram owicz et al . 1988) to describe the structure of the 
disk in both regimes. 

To calculate the disk's properties, we solve the equa- 
tions of conservation of mass, momentum, and energy: 



M = —AnRHpVr , 
3 v 1 



q =q 



pTVr 



R ' 



(15) 
(16) 

(17) 



where we have approximated the radial entropy gradient as 
ds/dR ~ —s/R. Here M is the accretion rate, R is the cylin- 
drical distance from the BH, H is the disk scale height, p 
is the density, v r is the radial velocity, and T is the mid- 
plane temperature. The no-torque boundary condition at 
the inner edge of the disk implies / = 1 — (-RlscV-R) 1 ^ 2 - 
We neglect gas pressure, since radiation pressure is domi- 
nant throughout the disk for at least a few years; we further 
assume that the visco us stress is proportional to the radia- 
tion pressure, so that (|Shakura fc Sunvaev| [l973) v — ac B H 
with sound speed c s = (aT 4 /ip) 1 ^ 2 and H — c s /Qk, where 
Qk = (GMbh/R 3 ) 1 ^ 2 - Simulations indicate that this as- 
sumpti on is reasonable a nd that such disks are thermally 
stable (|Hirose et al. 2009). The vertically integrated heating 
and cooling rates are given by 2Hq + = 3GA/bh Mf/itrR 3 
and 2Hq~ = 8crT 4 /3r, where the half-height optical depth 
is t — HpK s and k s is the electron scattering opacity. These 
relations form a quadratic equation for the dimensionless 
quantity K s aT 4 /cflK = 3t(c b /c), 







n B aT 4 



4 

"3a 



K s aT 4 
cQk 



1L 

3a 2 



10M 



(!) 



.(18) 



Solving equation (|18p yields the effective temperature of the 
disk, 



, 4aT 4 3GM BH M/ 

Crieff = — = ^~=^ X 



3r 



8ttR 3 



1 I 1 3 / WM 

2 + i 4 + 2 ; UfEdd 



(!) 



1/2 



(19) 



Combining this relation with equation ([2]), we calculate the 
luminosity and spectrum of the disk as a function of time, 
modeling it as a multicolor blackbody. 

The solution to (|18p also yields the disk scale height 
ratio, 



H_ 

~R. 



-J 



10A/ 



+ 



(.Rs) 



10M 

MEdd 



(r S ) 



1/2 



(20) 



The scale height H ~ R while Mf a iib ac k ^ A^Edd, and de- 
creases as t~ 5 ^ 3 at fixed R at later times. The viscous time 
at a radius R in the disk is 

-1/2 



fv 



GM B 



V R 3 



( 



(21) 




-0.08 



Figure 1. Spatial diagram of the equatorial stellar debris in the 
plane defined by the original orbit of the star. The unbound ma- 
terial is shown 0.1 yr, 0.3 yr, and 1 yr after the tidal disruption of 
a solar-type star by a 10 6 Mq BH at R p = 3Rs- The dashed and 
dotted curves indicate the incoming parabolic trajectory of the 
star and its continuation if the star were not disrupted. The de- 
bris also has an inclination dispersion perpendicular to this plane 
of Ai ~ 2_R*/_R P . 



which is ~ a" 1 times the local dynamical time during the 
super-Eddington phase, and later increases as t 10 ^ 3 at fixed 
R. For a ~ 0.1, the viscous timescale evaluated at the disk's 
outer edge is shorter than the time t since disruption for 
~ 1 — 3 yr; our assumption of steady-state accretion during 
this period is thus consistent. 



3 THE EQUATORIAL UNBOUND MATERIAL 

While half of the initial star becomes bound to the BH dur- 
ing the disruption, the other half gains energy and escapes 
from the BH on hyperbolic trajectories. From the viewpoint 
of the BH, this unbound material subtends a solid angle 
AQ, with a dispersion A<f) in azimuth and a dispersion Ai 
in orbital inclination. This material absorbs and re-radiates 
a fraction of the luminosity from the accretion disk0 We 
now estimate the dimensions of the unbound wedge. 

In the orbital plane at a fixed time t ;> tfaiibackj the 
unbound stellar debris lies along an arc, as the spread in 
specific energy produces a spread in radius and azimuthal 
angle (see Fig. [TJ . The most energetic particles escape on a 
hyperbolic orbit with eccentricity e max ~ l + GR+fRp. These 
particles race away from the BH at a substantial fraction of 
the speed of light, 



2 The polar outflow could also irradiate the unbound material, 
but it will have less of an effect because its spectrum is softer and 
its luminosity declines more rapidly. 
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R P 



1/2 



Q.3M, 



-1/2 



R 



p,3R s 



(22) 



(ignoring relativistic effects) and lie furthest from the BH at 
a distance 



Rn 



0.01M, 



-1/2 



-«p,3R s r * 



i 



0.1 yr 



pc . 



(23) 



They also have the smallest angle away from stellar pericen- 
ter, min ~ f x , where /oo obeys cosf x = -l/e max so tha t 
0mm ~ 7T-(12i?*/i? p ) 1/2 (see also lKhokhlov fc Melialfl996h . 
Particles with lower energies and thus smaller eccentricities 
are closer to the BH and make a larger angle relative to peri- 
center, up to 4> ~ This produces an azimuthal dispersion 
A<j) ~ {12R-,/R p ) 1/2 . 

Particles having the maximum vertical velocity A« p 
have vertical specific angular momentum j z = R p v p , and 
total specific angular momentum 



(24) 



to lowest order in Av p /v p . The orbital inclination i is given 



by cosi = ~ 1 - (l/2)(7?,/ J R p ) 



±R*/Rp. The 



resulting inclination dispersion is At ~ 2R ic /R p (our result 
is consistent with Evans & Kochanek 1989 but we disagree 
with Khokhlov & Melia 1996). 

The finite inclination dispersion produces a vertical wall 
of debris whose inside face scatters, absorbs, and re-radiates 
a fraction of the disk's emission. This face subtends a solid 
angle 



An = AiA0 



3/2 



48 V2 ( «* 
Ho 



0.2M- 3/2 R^gr 3 J 2 sr. 



(25) 



The number density of particles in the unbound wedge 
is n ~ {M i ,/2rn p )/(R 2 ARAtt/3), where Ai? is the radial 
dispersion of the material at fixed 0. This dispersion is 
due to differences in the particles' radial velocities and az- 
imuthal positions when the star passes through pericenter. 
Particles at ~ i? max travel on orbits whose pericenter is 
shifted from the star's pericenter by an angle ~ ±3(i?*/7? p ), 
which produces a spread in radial position (A_R/i?) max ~ 
(3-R* / Rp) 1 / 2 . The number density is then 



n ~ 10 9 Ml /2 Rl 



p,3R s m * r - 



-7/2 



0.1 yr 



cm 



(26) 



and the radial column density seen by the black hole is 

-2 



N ~ 10 25 M 5 6 /2 Rl /2 Rs m+r 



-5/2 



t 



0.1 yr 



(27) 



As the unbound material expands, it cools very quickly; 
after at most a few weeks, the gas would all be neutral if 
not for the disk's ionizing radiation. This radiation ionizes 
the surface layer of the unbound material. The ionized gas 
in turn emits via bremsstrahlung, radiative recombination, 
and lines. The physical conditions and processes here are 
similar to those in the broad line region of an active galactic 
nucleus (AGN). 

The ionized gas can reach photoionization equilibrium 
provided conditions change more slowly than the hydrogen 



recombination rate t'J, ~ n a a Icc . The recombination coeffi- 
cient for hydrogen is a ICC ps 4 x 10~ 13 cm 3 s -1 , and n c is the 
electron number density. In the ionized region, n c /n ps 1, as 
we show below. The material can remain in equilibrium for 
at least a few years, until t ICC /t ;> 1: 



tree , jA-l m-3 »#-7/2 n-5 -1 7/2 

— ~ (na IC ct) ~ 10 M 6 R 3R ra* ry 



t 



0.1 yr 



.(28) 



The column depth of the ionization front is Ai on ~ 
10 2i U cm -2 , where U = Ldi s k/47rJ? 2 nax c(/i^)n is the ioniza- 
tion parameter, 



U 



0.3 



(29) 



\L E dd) [ 0.1 



(hu) 



keV 



^ 6 " 3/2 ^p" 3 3R s ^ 1 ^ /2 



0.1 yr 



The electron density in the ionized layer is n e m n(l — 
10 _6 ?7 _1 ) ps n and the fractional depth of the ionization 
front is 
AR ion N ion 



AR N 

-3 f L disk \ ( (hv) 



3x10' 



/ ^diBk \ / JJ 

\L Edd J I 0.1 



M7^R. 



13/2 -2 5 



so the ionized layer is typically thin and highly ionized. 
4 PREDICTED EMISSION 

We use the results of §!j2] and [3] to calculate the emission 
due to the tidal disruption of a solar-type star as a func- 
tion of time and wavelength. We consider a solar-type star 
because the stellar mass-radius relation and typical stellar 
mass functions imply that these stars should dominate the 
event rate. The two key parameters we vary are the star's 
pericenter distance R p and the BH mass Mbh- We consider 
the mass range M B h ~ 10 5 - 10 s M . 

4.1 Super-Eddington Outflows 

Early on, when Mf a ii ba ck ^ MEdd (t < t E dd), outflowing gas 
likely dominates the emission. We calculate its properties 
using results from £12.11 

In Figure [2] we plot the spectral energy distribution at 
various times during the outflow phase, for three fiducial 
models: M BH = 10 6 Af Q and R p = 3R S ; Mbh = 10 6 M© 



t 



and _R D 



Rt; and Mbh 



1O 7 M and R p = R T . We 



take nominal values of /„ = 1 and / ou t = 0.1; we discuss 
the uncertainties in these parameters in ^5.31 and [6] The 
photosphere lies well inside the edge of the outflow at all 
times shown in Figure [5] The emission from the outflow 
has a blackbody spectrum, initially peaking at optical/UV 
wavelengths. As time passes and the density of the outflow 
subsides, the photosphere recedes and the emission becomes 
hotter but less luminous. 

In Figure [3] we plot g- (4770A) and i- (7625A) band 
light curves for the three fiducial models. In the leftmost 
panel at t <; 1 day, the edge of the outflow limits the size of 
the photosphere, so the photosphere initially expands, fol- 
lowing the edge of the outflow. After a time t e d ge (eq. [S]), 
however, the photosphere begins to recede inside the edge 



,(30) 



© 2009 RAS, MNRAS 00Q.[TlfT7l 



6 L. E. Strubbe and E. Quataert 




Figure 2. Spectral energy distributions for our three fiducial tidal disruption flares at several different times after disruption. These 
spectra include only the emission from the super-Eddington outflows (for / ou t = 0.1 and f v = 1; see eqns [5] & [6]), which dominate the 
emission at early times (see Fig. [3}. For Afgjj = 10 7 Mq and R p = i?x (right panel), the disk dominates the emission for t J> 50 days 
(Fig. [3j) ; this is why we do not plot the outflow emission at later times. 




Figure 3. Light curves at g- (solid) and i-band (dashed) for our three fiducial models. In each panel, times to the left of the vertical 
dotted line (~ *Edd; eq. [3]) have super-Eddington fallback rates and an outflow dominates the optical emission; at later times, the 
fallback rate is sub-Eddington and the emission is produced by the accretion disk and the photoionized unbound material (see £|4.2| |. 
In the leftmost panel, the emission rises at early times because the edge of the outflow limits the size of the photosphere. The optical 
emission then declines until the end of the outflow phase as the photosphere recedes and the photospheric temperature rises (Fig. [2). 



of the outflow and the luminosity declines. In the middle 
and rightmost panels, the photosphere lies well inside i? c d g c 
for virtually the entire outflow phase. The optical emission 
decreases as the photosphere's emitting area decreases and 
the temperature rises only slowly. As Figures [2] & demon- 
strate, the peak optical luminosity of the outflow is sub- 
stantial, ~ 10 43 ergs -1 ~ lO 9 Z/0, comparable to the optical 
luminosity of a supernova. The color of the emission is quite 
blue (g — r ~ —0.8). To illustrate how the peak luminosity 
depends on the parameters of the disruption event, Figure 
2] shows the peak bolometric and g-band luminosities of the 
outflow as a function of Mbh, for R p = 3Rs and R p = Rt- 



For sources at cosmological distances (which are detectable; 
fJS|, the negative k-correction associated with the Rayleigh- 
Jeans tail implies that the rest-frame g-band luminosity in 
Figure [4] underestimates the peak optical luminosity visible 
at Earth. 



4.2 Disk and Photoionized Unbound Debris 

When Mfaiiback i> Mvdd, a fraction of the falling-back gas is 
blown away while the remainder likely accretes via an advec- 
tive disk ( H2.2p . As the fallback rate declines below Edding- 
ton, the photons are able to diffuse out of the region close 
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M BH (M ) 

Figure 4. Peak bolometric (heavy lines) and (/-band (light lines) 
luminosities of the early-time super-Eddington outflows as func- 
tions of Mbh, for Ftp = 3-Rg (dashed) and R p = Rx (solid), 
assuming that / ou t = 0.1 and f v = 1. Figure [TT] shows the dura- 
tion of this phase. 



to the BH and the disk cools efficiently, but aiso becomes 
fess luminous. The vertical dotted line in Figure[3]delineates 
the super-Eddington fallback (and outflow) phase from the 
sub-Eddington fallback phase. 

The accretion disk irradiates the surface of the equato- 
rial unbound stellar material (331 . In this section we cal- 
culate the combined emission produced by the accretion 
disk and the irradiated stellar debris. In order to isolate 
the more theoretically secure emission by the disk and pho- 
toionized material, we do not consider the emission from 
super-Eddington outflows in this section. We show results 
for the disk and photoionized material at both t < t^dd and 
t > tEda; depending on the geometry of the outflow, and 
the viewing angle of the observer to the source, it is possible 
that all three emission components could be visible at early 
times. Because the mass driven away by outflows during the 
super-Eddington phase can also be photoionized by the cen- 
tral source at times t > t^dd, our emission line predictions 
are likely a lower limit to the total emission line fluxes (^6]). 

We calculate the photoionization properties of the un- 
bound material using version 07.02.0 2 of the public l y avai l- 
able code Cloudy, last described by iFerland et al.l l|l99Sl ). 
We simplify the geometry: the unbound spray traces out a 
widening spiral shape with most of the area coming from 
close to ~ -Rmax, so we approximate it as a cloud of area 
J? max AS7 located a distance i? max from the ionizing source. 
Our model cloud has constant density n (eq. [26]), col- 
umn depth iV (eq. [27]), and is irradiated by the accre- 
tion disk having the luminosity and spectrum described in 
£|2.2I The total emission calculated here is the sum of the 
emission from this photoionized layer and the emission from 
the central accretion disk. We focus on non-rotating BHs 



(Rlso = 3i?s), although we quote results for rapidly rotat- 
ing holes (J?lso = Rs) as well. 

In Figure [5] we plot the spectral energy distribution 30 
days, 100 days, 300 days, and 1000 days after disruption, 
for our three fiducial models: Mbh = 10 6 Mq and R p — 
3R S ; Mbh = 1O 6 M and R p = R T ; and Mbh = 1O 7 M 
and R p = Rt- The early-time short- wavelength peaks at 
~ 0.1 keV with luminosity ~ L^dd are emission from the 
disk. After a time tEdd, the mass fallback rate declines below 
the Eddington rate, and the disk begins to cool and fade. 
For R p ~ R LS o and Mbh ~ 10 5 - 1O 6 M , the optical light 
is dominated by lines and continuum from the photoionized 
material. For larger Mbh (and larger R p /Rs), the equatorial 
debris subtends a smaller solid angle (see eq. [35]) and the 
disk's luminosity is larger, so the disk dominates the optical 
emission. 

Figure [5] zooms in on the UV/optical/near-infrared 
spectra for our three fiducial models. The emission lines 
characteristic of the broad line region of an AGN are 
typically the strongest lines here as well: e.g., Lya, CIV 
1548+1551, H/3, and Ha. In most cases, these lines are 
optically thick for more than a year. The lines are ex- 
tremely broad, since the marginally bound gas has a speed 
close to zero while the most energetic gas leaves the BH at 
•Umax ~ 0.4c, 0.09c, and 0.2c for Mbh = 1O 6 M , R p = 3R S ; 
Mbh = 10 s M , R p = R T ; and Mbh = 1O 7 M , R p = R T , 
respectively (see eq. [22]). In addition, the mean velocity 
along our line of sight will usually be substantial, so the 
lines should have a large redshift or blueshift on top of the 
galaxy's redshift. For clarity, we plot the spectra with a mean 
redshift of zero. 

Figure [7] focuses on the evolution of five strong lines, 
plotting the ratio L^n^/L^coat at line center for each0 The 
quantity L^iine is the line intensity at line center account- 
ing for the significant broadening. As the surface area of the 
equatorial wedge grows in time, line luminosities grow un- 
til Mfaiiback <; M e dd and irradiation by the disk subsides. 
The quantity Z/„, C ont is the continuum intensity, which in- 
cludes the contributions of both the disk (blackbody) and 
the photoionized unbound material (bremsstrahlung and ra- 
diative recombination) — again, the emission from the super- 
Eddington outflows is not included in L v , con t. The lines re- 
main prominent for a few years, and are strongest and broad- 
est for small Mbh and small R p /Rs. The UV lines are the 
strongest lines when the unbound material dominates the 
continuum (left panel), while the near-infrared lines are the 
strongest when the disk dominates the continuum (middle 
and right panels). 

We next describe the broadband optical evolution of 
a tidal disruption event. Figure [3] also discussed in S|4. 1 1 
plots the optical light curve for each fiducial model, show- 
ing the total emission at both g- and i-bands. For t <J t^dd 
(left of the dotted lines), the emission is dominated by the 
super-Eddington outflows, while for t ;> t^dd the emission is 
dominated by the accretion disk and photoionized equato- 
rial debris. Once t > t E dd, i.e., Mfaiiback <; Medd, the disk's 
optical luminosity falls off gently, approximately as t~ 5 ^ 12 : 



3 By comparing the results of different versions of Cloudy, we 
find that the results for line strengths can be uncertain by up to 
factor few. 
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Figure 5. Spectral energy distributions for tidal flares around a non-rotating BH, 30, 100, 300, and 1000 days after disruption. Emission 
from the accretion disk dominates at short wavelengths. The photoionized unbound stellar debris absorbs and re-radiates some of the 
disk's emission, producing optical-infrared emission. These spectra do not include the emission from super-Eddington outflows at early 
times; see Figs. [2] & [3] for this emission. 




1000 10000 1000 10000 1000 10000 

A (Angstroms) A (Angstroms) A (Angstroms) 

Figure 6. UV to near-infrared spectra for tidal flares around a non-rotating BH, 30, 100, 300, and 1000 days after disruption. The 
spectra are the sum of contributions from the accretion disk and the photoionized unbound material, but do not include the emission 
from super-Eddington outflows, which likely dominate at early times (Fig. |3J. The linewidths and line strengths are both larger for 
smaller R p /Rs and smaller A/bh (eq. [22]). 



although the bolometric luminosity is declining as t~ ' , the 
optical emission lies on the Rayleigh- Jeans tail. Increasing 
Rp/Rs and/or Mbh increases the disk's luminosity by up to 
two orders of magnitude because of the disk's larger emit- 
ting area and/or because LEdd rises. At all times, the disk 
emission is quite blue (g — r ~ — 1). 

For large Mbh and/or large R p /Rs (middle and right 
panels in Fig. [3}, the disk outshines the photoionized ma- 
terial at optical wavelengths, and the light curves and color 
evolution are determined by the disk emission alone. By con- 
trast, for M BH = 1O 6 M and R p = 3R S (left panel), the 
photoionized material's optical line emission is initially an 
order of magnitude brighter than the disk. As the illumi- 



nating power of the disk declines but the unbound debris 
becomes less dense, different lines wax and wane. The sig- 
nificant redshift or blueshift of the unbound material fur- 
ther complicates the photometry by altering which lines 
contribute in which wavebands (again, our figures assume 
a mean redshift of zero). These effects can produce a non- 
monotonic light curve and a complicated color evolution, 
depending on the exact redshift of the source and the veloc- 
ity of the equatorial debris. 
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Figure 7. Line strength evolution for the three fiducial models. Results are shown for Lyc* (solid), CIV 1548+1551 (dotted), H/3 (dashed), 
Ha (dot-dashed), and Pact (triple-dot-dashed). The quantity y lne / L u>con t is the ratio at line center of the line intensity to continuum 
intensity, taking into account emission from both the disk and the photoionized material, but not the early-time super-Eddington outflows. 



Table 1. Assumed parameters for transient surveysSl and predicted rates. Our results can be scaled to other survey parameters and 
model assumptions using eq. I|32|l . 



Survey 


/lim 


/sky 


Cadence 


Rate: D+UB a (yr" 1 ) 


Rate: Outflows'' (yr _1 ) 


Pan-STARRS 3tt Survey 


23 AB mag (g, i-band) 


0.75 


6 months 


4-12 


200 


Pan-STARRS Medium Deep Survey (MDS) 


25 AB mag (cj-band) 


2 x 10" 3 


4 days 


0.2 - 1 


20 


Palomar Transient Factory (PTF) 


21 AB mag (g-band) 


0.2 


5 days 


0.3 - 0.8 


300 


Large Synoptic Survey Telescope (LSST) 


24.5 AB mag (g-band) 


0.5 


3 days 


60 - 250 


6000 


Synoptic All-Sky Infrared Survey (SASIR) C 


23.5 AB mag (Y-band) 


0.03? 


10 days? 


0.1 - 0.5 


100 


ROSAT All-Sky Survey d 


2 X 10~ 12 ergs" 1 cm" 2 


1 


e 


2 - 100 


N/Af 


GALEX Deep Imaging Survey 


25 AB mag (2316 A) 


7 x 10" 5 


e 


0.05 - 0.2 


N/A/ 



a Rates for the emission from the disk and unbound equatorial debris; the range corresponds to i?LSO = (l° w )i ^LSO = 1-Rs 
(high). This emission is relatively faint in the optical/infrared and may be difficult to detect relative to the host bulge (see ilBl . 
b Rates for the emission from super-Eddington outflows, restricted to z < 1, for / ou t = 0.1 and f v = 1 (see Fig. 1 1 3 D . 

c The survey strategy for SASIR has not yet been finalizerL 

d These parameters are for comparison with the all-sky rate calculation by iDonlev et al. | d2002h . 
e ROSAT and GALEX do not have regular cadences. 
f ROSAT and GALEX have insufficient cadence and/or sky coverage to detect flares from super-Eddington outflows. 



5 PREDICTED RATES 

We use our calculated spectra and light curves to predict 
the number of tidal disruption events detectable by obser- 
vational surveys. We focus on an (almost) all-sky optical 
survey like the Pan-STARRS PS1 3n survey, but we also 
predict results for surveys with more rapid cadence (e.g., 
PTF and LSST) and discuss the results of our models com- 
pared to ROSAT and GALEX observations. Our assumed 
survey parameters are listed in Table [T] along with some 
of our results Our results can readily be scaled to other 
surveys using equation (|32[) discussed below. 

To predict rates, we use the red s hift-de pendent BH 
mass function given by iHopkins et all i|2007f ). At z ~ 0, 



4 These are intended to be illustrative, and may not correspond 
precisely to the true observational survey parameters, although 
we have attempted to be as accurate as possible. 



the BH density is ~ 10" 2 Mpc" 3 for - 10 6 M Q and gently 
falling at higher masses; as z rises to 3, the BH number den- 
sity falls by ~ 30. We assume that the BH number density 
for 10 5 - 1O 6 M is the same as for ~ 10 6 M Q , although it is 
poorly constrained observationally. We do not consider tidal 
disruption events beyond z ~ 3. (Our results can easily be 
scaled to other assumed BH mass densities; see eq. [32]). 



The rate of tidal disruptions within a single galaxy 
is 7(A/bh). To predict detection rates, we assume that 7 
is indepe ndent of BH mass. W e adopt 7 = 10~ 5 yr" 1 as 
found by IDonlev et ail (|2002h using the ROSAT All-Sky 
Survey, which is also in line with conservative theoretical 
estimates. We further assume that this rate is distributed 
equally among logarithmic bins of stellar pericenter distance 
-Rp, so that d7/dln7ip = 7/ ln(ii/r/-Rp,min). In the limit of 
z < 1, the equation for the predicted rate is 
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dT 



din Mesh 



4tt 



d m ax/sky 



dn 



(h. 



din A/bh din R p 



dlnR p (31) 



where / s k y is the fraction of sky surveyed; when necessary we 
use the generalization of equation (|31[1 that includes cosmo- 
logical effects. When the duration of a flare tn aTC is shorter 
than the cadence of the survey t ca d, we approximate the 
probability of detection as tflare/icad- 

We start by considering emission from only the ac- 
cretion disk and photoionized equatorial debris. Then in 
i !5.3l we include the emission from super-Eddington out- 
flows, where the physics is somewhat less certain, but the 
observational prospects are particularly promising. 



5.1 Disk and Photoionized Material 

For all but the largest Mbh, the duration of peak optical 
emission for the accretion disk and photoionized material is 
iflarc ~ ^Edd (eq. [3]). This timescale depends on the BH's 
mass and the star's pericenter distance, as shown in Fig- 
ure [8] For R p ~ Rt, the flare lasts for tEdd ~ 0.3 — lyr 
and then decays only gently since the disk dominates the 
optical emission. However, for M BH ~ 10 5 - 10 6 A/ Q and 
Rp ~ Rlso, the optical flare is shorter, tEdd ~ 0.03 — 0.1 yr, 
and then the emission decays more quickly since irradiation 
of the unbound material — the dominant source of optical 
emission — subsides. For Mbh ^8x 10 7 Mq, ta^ c ~ ifaiiback 
(eq. PP) because the fallback rate in these systems is never 
super-Eddington. 

Figure [9] shows our calculated rates for optically- 
detected tidal flares (for a survey like the Pan-STARRS 37r 
survey), for both non-rotating and rapidly rotating BHs. 
For Mbh <; 10 6 M , the disk contributes most of the emit- 
ted power, so the rates increase with Mbh as Ldisk ~ ^Edd 
increases. The rates are dominated by R p ~ Rt and Mbh ~ 
2 x 1O 7 M (non-rotating BHs) and Mbh ~ 10 8 M Q (rapidly 
rotating BHs). Since most of the flares that dominate the 
rates have relatively long durations (Fig. [8}, imperfect sur- 
vey cadence only modifies the detection rates by ~ 50%. At 
10 5 — 1O 6 M0 and small R p , the photoionized material re- 
emits a relatively large fraction of the disk's power in the 
optical and boosts the detection rates significantly. 

Integrated over A/bh = 10 6 - 10 s M , our estimated 
rates for the Pan-STARRS 3n survey, assuming non-rotating 
BHs, are 4yr _1 and 2yr _1 in g-band and i-band, respec- 
tively. The mass range 10 5 — 10 6 Mq contributes another 
0.4 yr -1 (both g-band and i-band), assuming (probably op- 
timistically) that dn/d In Mbh and 7 are the same at 1O 5 M 
as at 1O 6 M . If the BH is rotating faster, R p /Rs can be 
smaller. This allows an accretion disk to form for even 
A/bh ~ few x 10 7 - 10 s M , widens the disk for all Mbh, and 
increases the solid angle of the unbound material. Indeed for 
Mbh ~ 1O 5 M and R p ~ Rs, the unbound material covers 
a quarter of the sky! (At this point Rs ~ Rq so our approx- 
imations begin to break down.) These effects raise the total 
predicted rates for rapidly rotating BHs to ~ 12 yr" 1 . 

Figure [10] plots the detection rates as a function of 
R p /Rs for Mbh = 10 5 A/ o , lO 6 A/ , and 1O 7 M , for the 
disk alone (light lines) and for the disk plus photoionized 
material (heavy lines). The photoionized material enhances 
detection rates significantly for most R p /Rs at 10 M and 
1O 6 M , but has little effect for 1O 7 M . The rates decrease 



1000 




M BH (M ) 



Figure 8. Duration of maximum luminosity during the late-time 
accretion disk phase, plotted as a function of A/bh f° r Rp = 3/?s 



and R„ 



Rrp. Most flares last longer than ~ a month. The 



super-Eddington outflow at early times produces a shorter flare 
that should precede this emission (Figs.l3lfe llHfT3t . (The fallback 
rates for Mbh J> 8 X 1O 7 M are never super-Eddington, so the 
flare duration in these systems is ~ ifaiiback-) 



substantially for R p /Rs — > 1.5 because then the outer ra- 
dius of the disk ~ 2R P — 3Rs is at the last stable orbit; 
our disk model assumes a no-torque boundary condition at 
7?lso, implying that there is essentially no emission from 
the disk when R p ~ -Rlso- 

Although the rates quoted above and in Figures [9] and 
[I0]are for a survey covering 3/4 of the sky, assuming constant 
7(A/bh) and constant BH mass density below ~ 1O 6 M , our 
predicted rates can be scaled to other assumed parameters: 



dT 



din A/bh 



^ /lim /sky /cad 



dn 



d In A/bh 



7 (A/bh) 



(32) 



assuming that the sources are at z <C 1, where / ca d = 
min(f flarc /t cad , 1). 

The Pan-STARRS MDS, PTF, and LSST will have ca- 
dences of only a few days, giving them sensitivity to all 
events within their survey volumes. However, for detect- 
ing this relatively long-lived emission, the advantage of fast 
cadence is only minor; more significant is the spatial vol- 
ume probed by the survey (see Table [IJ. The MDS will 
deeply (biab ~ 25) image a relatively small (84-deg 2 ) re- 
gion, and so will detect ~ 5 — 10% as many events as the 
37r Survey, ~ 0.2 — 1 per year. PTF will image a larger 
region (8000 deg 2 ) but less deeply (mAB ~ 21), and so 
will have detection rates similar to the MDS. The emis- 
sion from early-time super-Eddington outflows may signifi- 
cantly increase these rates, as we discuss in £15.31 In the next 
decade, LSST will image a large region (20,000 deg 2 ) deeply 
(mAB ~ 24.5), and so will discover hundreds of tidal flares. 
MDS and LSST will detect events at cosmological distances 
(z m ax ~ 0.3 — 0.6), where the negative k-correction of the 
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Figure 9. Predicted detection rates as a function of BH mass for the Pan-STARRS 37r survey (see Table [TJ. Results are shown for 
non-rotating BHs (left panel) and rapidly rotating BHs (right panel). The rates shown here can be scaled to other surveys and other 
model parameters using eq. (132 D . These rates do not include the emission from super-Eddington outflows at early times; see Figs. Il2l fc 
I13l for these results. 



disk's blackbody peak enhances rates by a factor of a few. 
By co-adding images up to ~ 1 month, these short-cadence 
surveys will also be able to raise their detection rates of 
events having Mbh ~ 10 6 — 10 s Mq by a factor of a few. Also 
in the next decade, SASIR plans to deeply image 140 deg 2 
each night in the near infrared (to mAB ~ 23.5 at Y-band; 
iBloom et al.ll2009h . Their observing strategy is not yet final- 
ized; assuming a 10-day cadence covering 1400 deg 2 (a good 
strategy for detecting flares during the super-Eddington out- 
flow phase; see ^5.3f) . this survey should detect a flare every 
few years from the accretion disk plus photoionized debris. 



5.2 Current Observational Constraints 

As a check on our model, we have calculated rates for the 
ROSA T All-Sky Survey for comparison with iDonlev et al.l 
(2002 )'s result of 42 events per year over the whole sky for 
/i im = 2 x 10" 12 ergs -1 cm -2 (0.2 - 2.4 keV)Q Considering 
blackbody emission alone, we predict 2yr _1 (R\ B = 3-Rs) 
or 40yr _1 (Ri n = Rs)- This strong sensitivity to R ln arises 
because the ROSAT band is on the Wien tail of the disk 
emission. If we assume that 10% of the emission is in an 
X-ray power-law tail with a photon index V = 3 (not un- 
reasonable as^rnirrjtions_Jor XjTay emission from accreting 
BHs; e.g.. iKoratkar fc Blaeslll999h . our predicted rates are 
10 yr" 1 (R ia = 3R S ) and 100 yr" 1 (R in = R s ). Given the 
large uncertainties in the X-ray emission, our predictions are 
consistent with the observational results. We also compare 
with detection rates in the GALEX Deep Imaging Survey 



5 The cadence is ~ 1 year but irregular; we assume it is perfect 
for simplicity. The rates we predict thus may be slightly high. 



l|Gezari et al.ll2008h . They search an area of 2.882 square de- 
grees, observed at FUV (1539 A) and NUV (2316 A) down 
to /ii m ~ 25 AB magnitudes. Gezari et al. detect 2 eventijf] 
in this region over 3 years. This is somewhat higher than 
our predicted rates in the NUV, 0.05 yr _1 if Ri n = 3Rs and 
0.2yr _1 if _Ri n = Rs (assuming perfect cadence for simplic- 
ity), and may suggest that the disruption rate per galaxy 7 
is a factor of a few higher than we have assumed here. 

As the above estimates demonstrate, consistency with 
GALEX and ROSAT constraints prefers a rate per galaxy 
of 7 ~ (1 — 3) x 10 _5 yr _1 . Significantly l arger disruption 
rates , as some calculations predict (e.g., iMerritt fc Poonl 
2004), are inconsistent with current observational limits un- 
less dust obscuration has significantly biased the ROSAT 
and GALEX results or the large disruption rates are con- 
fined to brief epochs in a galaxy's life (e.g., a merger). Note 
also that this constraint only applies to massive BHs with 
M BH ~ 10 7 M o - 10 8 M Q , because UV and X-ray surveys 
select for these systems ( !j6j> - 

5.3 Super-Eddington Outflows 

We now calculate optical detection rates accounting 
for emission from the luminous but shorter-lived super- 
Eddington outflows. We begin by plotting the duration of 
peak emission for the outflow phase, tfl ar o ~ ^outflow, in Fig- 
ure[TT]as a function of Mbh for R p = 3Rs and R p — Rt- For 
R p — 3Rs and Mbh ^6x 10 6 M Q , the duration of the out- 
flow is set by the time at which the photosphere recedes in- 
side the edge of the outflow (t ou tflow ~ iedge! eq. 0) while for 
larger A/bh, the duration is set by the timescale for the most 

6 The candidate of lGezari et al. I j2009l) is in a different field. 
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Rp/Rs 

Figure 10. Predicted detection rates as a function of pericenter 
distance R p /Rg for the Pan-STARRS 3tt survey at g-band (see 
Table [TJ. Results are shown for M BH = 1O 5 M (solid), 10 6 M Q 
(dashed), and 1O 7 M (dot-dashed), assuming a non-rotating BH. 
The thin/light lines are for the disk emission alone, while the 
thick/heavy lines include the emission from both the disk and 
the photoionized unbound material. The photoionized material 
significantly increases the rates for low Mbh and small R p /Rs- 
These rates do not include emission from supcr-Eddington out- 
flows at early times; those results are shown in Figs. [T2lfc 1131 The 
small fluctuations in the results for Mbh = 10 5 Mq are due to dif- 
ficulties in performing Cloudy calculations at very high densities. 




10 5 10 6 10 7 10 8 
M BH (M ) 

Figure 11. Duration of peak luminosity during the early super- 
Eddington outflow phase, as a function of Mbh f° r Rp = 3i?g 
and R p = i?x (f° r /out = 0.1 and /„ = 1; see eqs. [5] & [6]). The 
vertical dotted lines mark the boundary between events where 
the edge of the outflow limits the size of the photosphere (lower 
BH masses) and where it does not (higher BH masses). The flares 
from super-Eddington outflows typically last longer than the few- 
day cadence of surveys like PTF, Pan-STARRS MDS, and LSST, 
but they are often short enough that they would not be detected 
in surveys optimized solely for supernovae. 



bound material to return to pericenter (t ou tflow ~ ifaiiback; 
eq. []]); for R p = Rt, the transition from t edge to tailback 
occurs at a somewhat lower BH mass of M BH £ 4x 1O 5 M . 
Figure [11] shows that most flares last longer than the few- 
day cadences of surveys like PTF, Pan-STARRS MDS, and 
LSST, but are much shorter than the months-long cadence 
of the Pan-STARRS 3tt survey. ROSAT and GALEX are 
unlikely to have detected events during this phase due to 
insufficient cadence and sky coverage. 

In Figure 1121 we plot detection rates as a function of 
Rp/Rs for M BH = 1O 5 M , 1O 6 M , and 1O 7 M , all for the 
Pan-STARRS 37r survey at g-band. At small R p /Rs, the 
edge of the outflow limits the size of the photosphere. As 
Rp/Rs increases, the radius of the photosphere at i e d g e in- 
creases and so the rates increase. For Mbh = 1O 6 M and 
1O 7 M at the largest R p /Rs, the photosphere is no longer 
limited by the edge of the outflow, and the maximum lumi- 
nosity occurs at ifailback- For large R p /Rs, the rate declines 
as the photosphere recedes inward. 

In Figure [131 we plot overall detection rates at g-band 
for Pan-STARRS 3tt and MDS, PTF, and LSST. In the left- 
most panel, we assume that 10% of the falling back material 
flows out in the wind (/ ou t = 0.1), as we have previously. 
Here we restrict detections to redshifts z < 1, where our as- 
sumed tidal disruption rate per galaxy is appropriate. The 
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Figure 12. Predicted detection rates for emission from super- 
Eddington outflows at early times (for / ou t = 0.1 and f v = 1; 
see eqs. [5] & [6]). Results are shown as a function of Rp/Rs 
for A/ BH = 1O 5 M , 1O 6 M , and lO 7 A/ , all for a Pan-STARRS 
37r-like survey at g-band (see Table [TJ- 
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Pan-STARRS 3tt survey should detect 200 yr -1 , while the 
MDS should detect 20 yr -1 because of its smaller spatial 
volume. PTF should detect 300 yr" 1 as well, since its fast 
cadence makes up for its smaller spatial volume relative to 
the Pan-STARRS 3n survey. LSST's large spatial volume 
and rapid cadence should allow it to detect 6000 yr~ ! As- 
suming the survey parameters and strategy described in E I5.1I 
and Table [1] SASIR should detect ~ 100 yr -1 as well. 

In the central panel of Figure [131 we plot detection rates 
for redshifts z > 1 to highlight the possibility of studying 
tidal disruption events at cosmological distances. Deep sur- 
veys like MDS will be sensitive to events far beyond z ~ 1. 
Even the shallower 37r survey could detect ~ 30 events per 
year at z > 1. Our predicted rates at z > 1 are particu- 
larly uncertain: higher redshift sources are more likely to 
be obscured, and the lower mass BHs of interest for tidal 
disruption are still growing s ignificantly in mass at these 
redshifts (|Heckman et al.ll2004T ) so the tidal disruption rate 
and mass function become less certain. Nonetheless, the fact 
that tidal disruption flares will be detectable at z > 1 with 
forthcoming surveys highlights that these sources may be- 
come a strong probe of the evolution of ~ 10 6 — 10 s Mq 
BHs. 

In the rightmost panel of Figure 1131 we show the sensi- 
tivity of our predictions to uncertainty in the outflow model 
by plotting the rates for the Pan-STARRS 37r survey for dif- 
ferent assumptions about the fraction of the material that is 
blown away in the outflow, /out - For / out ~ 0.3, the detection 
rate is 600 yr -1 , while for / Dut ~ 0.01 it falls to 8yr _1 . The 
rate falls rapidly at the highest Mbh because the outflow is 
optically thin. Figure [13] shows that even if we very conser- 
vatively assume that only 1% of the material is blown away 
when the fallback rate is super-Eddington, upcoming opti- 
cal surveys like Pan-STARRS should still be able to detect 
a significant number of tidal disruptions during the super- 
Eddington phase. 



6 DISCUSSION 

We have calculated the spectra and light curves produced 
by the tidal disruption of a solar-type star by a massive 
black hole. Upcoming optical transient surveys should de- 
tect many such events (35}. Our results demonstrate that 
there are at least three different emission components that 
are important during tidal flares: (1) outflows at early times 
when the fallback rate is super-Eddington, (2) a compact 
(<J 10 — 100i?s) accretion disk around the BH, and (3) stel- 
lar debris that is unbound during the disruption and forms 
an outflowing "wedge" in the equatorial plane (see Fig. [l}. 
It is also possible that the super-Eddington fallback powers 
a lower-density magnetically-dominated jet, but the proper- 
ties of such a jet are difficult to predict so we do not consider 
this potential source of emission. 

Each of these three components contributes to the to- 
tal emission from tidal disruption flares. At early times, the 
super-Eddington outflows likely dominate, producing a few- 
to 10-day optical-infrared flare with a luminosity compara- 
ble to that of a supernova fFigs. I2ll4l fc ITTjl. As the fallback 
rate decreases below the Eddington rate, these outflows will 
diminish, revealing the underlying accretion disk that emits 
primarily in the UV to soft X-rays (Fig. [5}; at this time, the 



optical emission is likely to be much less than that of a typi- 
cal AGN — and well below that of the super-Eddington phase 
(Fig. [3} — because the accretion disk is not very spatially ex- 
tended. The central UV and soft X-ray source photoionizes 
the inner edge of equatorial stellar debris, producing a spec- 
tum of broad emission lines (Figs. [5] [6] & whose "rest" 
wavelength should be either blueshifted or redshifted with 
respect to the host galaxy depending on the line of sight of 
the observer relative to the escaping material. We find that 
this spectroscopic signature of tidal flares is the strongest for 
low mass BHs because the equatorial stellar debris occupies 
the largest solid angle in these systems (eq. [25|1. 

Although the above stages are the focus of this paper, 
for completeness we briefly discuss the rest of the evolution 
of a tidal disruption event. As the fallback rate continues to 
decrease below the Eddington rate, the viscous time in the 
thin disk (eq. J2TJ) increases and becomes comparable to the 
time t since disruption; at this point, matter begins to build 
up rather than rapidly accreting onto the BH. For disrup- 
tion at R p ~ Rt, which is likely to produce a significant 
fraction of the events (Figs. [T0l fc ll2[) . we estimate that this 
occurs ~ 3 years after disruption, when the fallback rate 
has decreased to Mfaiiback ~ 0.2 ag 5 / 7 M 6 ~ 2 ^ 3 Msdd (where 
a = O.lao.i is the dimensionless viscosity). Separately, we 
expect a significant change in the thermodynamics of the 
disk when Mfaiiback — a 2 MEdd ~ 10" 2 MEdd- Below this 
accretion rate, the material will no longer cool efficiently 
when it circularizes and shocks upon returning to pericen- 
ter. Instead of cooling to form a thin disk, the material will 
be shock heated to form a geometrically thick, radiatively 
inefficient accretion flow. In general, both geometrically 
thin/optically thick and geometrically thick/radiatively in- 
efficient disks appear to be stable accretion solutions for a 
given M <J 10 _2 A/Edd- However, in this case, the boundary 
condition that matter shocks up to the virial temperature 
upon returning to pericenter picks out the radiatively in- 
efficient solution once Mfaiiback <, 10~ 2 A/Edd- Because the 
viscous time in a thick disk is ~ a -1 times the dynamical 
time, once the accretion becomes radiatively inefficient, the 
viscous time is always much shorter than the orbital period 
of matter returning to pericenter. Moreover, the transition 
to radiatively inefficient accretion happens at a time compa- 
rable to when matter would otherwise begin accumulating 
in a thin disk, particularly for more massive BHs. This sug- 
gests that there is typically only a limited range of accretion 
rates ( and time) during wh ich the "spreading disk" solu- 
tion of lCannizzo et all (|l990l ) applies. Instead, at late times 
matter will rapidly accrete via a thick disk and the accretion 
rate will decay as ~ Mfaiiback o c t~ 5 ^ 3 . As in X-ray binaries 
jRemillard fc McClintock1l2006l '). we expect that the thermo- 
dynamic transition at ~ 10 _2 MEdd will be accompanied by 
a significant change in the luminosity and spectrum of the 
disk, and perhaps also by the production of relativistic jets. 
This should be explored in more detail in future work. 

Having summarized our key results and the timeline of 
a tidal disruption event, we now discuss some uncertainties 
in our models, observational challenges to detecting tidal 
flares, and the astrophysical applications of studies of tidal 
disruption events. 
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Figure 13. Predicted detection rates for emission at g-band from super-Eddington outflows for various optical transient surveys. In the 
leftmost panel, we use our standard outflow model in which /out = 0.1 and f v = 1 (see eqns [5] & [6]), and restrict detections to z < 1. In 
the central panel, we plot detection rates for z > 1 to illustrate the potential for surveys to be sensitive to the cosmological evolution of 
BHs; note that PTF is not deep enough to detect flares beyond z ~ 1. In the rightmost panel, we vary / ou t from 0.3 to 0.01 to illustrate 
the sensitivity of our results to uncertainties in the outflow physics. 



6.1 Super-Eddington Outflows 

In fl2.ll we described our simple model for outflows driven 
when the fallback rate is super-Eddington. Energy conserva- 
tion implies that the falling back material initially remains 
bound after returning to pericenter and circularizing, but 
even small amounts of accretion can release additional en- 
ergy and drive a powerful outflow. There is, however, a sig- 
nificant uncertainty in precisely how much of the falling back 
material is blown away, and in the kinematics of the outflow. 
We assume that the gas expands roughly spherically from 
the BH, but the flow is probably somewhat collimated along 
the pole, due to the original angular momentum of the stel- 
lar debris. Some photons can then leak out through the sides 
of the outflow rather than continuing to drive the expansion; 
in this case, the overall emission would be somewhat hotter 
and fainter, with a dependence on viewing angle. 

We have parameterized the terminal velocity of the out- 
flow using Uwind = f v VoEc(2R p ), and the mass outflow rate 
using A/ ou t = /out Mf a iiback; in all of our calculations, we have 
assumed that /„ ~ 1. If the gas actually expands more slowly 
(f v < 1), its density will be larger, so the photosphere will 
be larger, increasing the optical fluxes and detection rates. 
In the extreme case in which there is no unbound outflow, 
but super-Eddington fallback leads to a radiation-pressure- 
supported atmosphere around the BH that slowly expands 
until the photons can diffuse out, we also expect significant 
optical luminosities during the super-Eddington phase (e.g., 
lLoeb fc Ulmerlll997r i. 

It is worth noting that our predictions for the radia- 
tion from super-Eddington outflows are particularly uncer- 
tain for low Mbh and small R p /Rs, when the edge of the 
outflow limits the radius of the photosphere and determines 
both the peak luminosity (Fig. [3J) and duration (Fig. 1 1 1 f) of 
the flare; this does not, however, significantly influence our 
total predicted rates (Fig. 1 12p . 

Figure [13] shows that even if the outflow rate is just 



a few percent of the fallback rate (/out ~ 0.01 — 0.1), the 
outflowing gas is sufficiently bright in the optical that forth- 
coming surveys should detect a significant number of tidal 
disruption flares. This makes early-time optical flares from 
the tidal disruption of stars an extremely promising candi- 
date for current and future optical/infrared transient sur- 
veys. (We discuss some practical issues associated with de- 
tecting these sources in ^6.30 In this context, we note that 
in the radiation hydrodynami c simulations of acc retion at 
M = 100Af E dd carried out bv lOhsuga et all l|2005l ). ~ 10% 
of the gas becomes unbound, ~ 10% accretes, and the re- 
maining ~ 80% is marginally bound and may (or may not) 
eventually accrete as well. These precise values will depend 
on the pericenter of the star, with smaller R p /Rs likely lead- 
ing to smaller /out, i.e., a smaller fraction of the gas being 
blown away. Future observational constraints on the lumi- 
nosity, spectrum, and timescale of the super-Eddington out- 
flow phase should be able to strongly constrain the value of 
/out in individual events. These results will have important 
implications for how massive BHs grow. In particular, if / ou t 
is typically modest, this would imply that black holes can 
accrete at rates far above the Eddington rate, perhaps help- 
ing to explain how supermassive BHs (A/bh Jj 10 8 Mq) can 
be observable as luminous quasars as early as z ~ 6. 

Figure [13] also demonstrates that deep optical surveys 
such as the Pan-STARRS MDS and LSST will be sensitive to 
tidal flares at high redshift. These surveys may thus provide 
a powerful probe of the BH mass function and stellar dynam- 
ics in galactic nuclei as a function of redshift. For example, at 
z ~ 0.1, BHs h aving Mbh <^ 10 7 Mq are still growing signifi- 
cantly in mass ( Hcckman ct al. 20(3) and thus their disrup- 
tion rates may evolve significantly with redshift. In addition, 
galaxy mergers, which are more common at z ~ 1 — 2, could 
substa ntially increase the tidal disruption rate: IChen et al.l 
(2009) find rates of up to ~ 1 yr _1 for ~ 10 5 years after the 
merger due to three-body interactions between stars and a 
binary BH. 
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6.2 The Accretion Disk, Photoionized Gas, & 
Broad Emission Lines 

We now consider several aspects of our model for the ac- 
cretion disk, broad emission lines, and the material un- 
bound during the disruption f§i |2.2l and [3)). Our accretion 
disk model is designed to describe the emission from the 
time when the disk first forms through the following few 
years. During much of this period, the fallback rate is super- 
Eddington, and we expect the disk to be optically and ge- 
ometrically thick, with radiation pressure dominating gas 
pressure. Once the fallback rate becomes sub-Eddington, 
the disk becomes geometrically thin and may be subject 
to viscous instab ilities (although it is thermally stable; 
iHirose et al. I I2009D . These instabilities may lead to addi- 
tional time dependence not captured in our models, par- 
ticularly at late times when the viscous time in the thin 
disk becomes comparable to the orbital period of the ma- 
terial falling back to pericenter. As described above, once 
Af fallback <J a 2 A/Edd, the density is sufficiently low that 
the flow becomes radiatively inefficient and our model is 
no longer appropriate. The disk will then heat up and its 
spectrum will become significantly harder. This phase may 
be detectable by hard X-ray transient surveys like EXIST. 

It is also unclear how much tidal forces will spin up the 
rotation of the star as it approaches pericenter. We have as- 
sumed that the star is maximally spun up, so that stellar de- 
bris is accelerated to relative velocities Av ~ v p (R i ,/ Rp) 1 ^ 2 
in the azimuthal direction. If in fact the spin-up is less ef- 
fective, the onset of the flare, which occurs at the fallback 
time (eq. [I]), will be later by a factor of few and the solid 
angle subtended by the unbound equatorial debris will be 
somewhat smaller. This will not change our qualitative con- 
clusions, only some of our quantitative results. 

It is important to stress that the line emission we pre- 
dict may well be an underestimate in all cases: the accretion 
disk will also photoionize the back edge of the material that 
was blown away during the super-Eddington phase, which is 
far from the BH once the outflows subside. Simple estimates 
indicate that the density and velocity of this outflowing gas 
are similar to that of the gas unbound at the time of disrup- 
tion; as a result, irradiation of this gas will produce addi- 
tional broad hydrogen lines. The equivalent width of these 
lines depends on the solid angle subtended by the super- 
Eddington winds, which, although uncertain, is likely to be 
significant. These lines are unlikely to depend as sensitively 
on Mbh and R p /Rs as the emission lines from the equatorial 
debris (see Fig. [6] for the latter). As a result, observations 
of the line emission will help constrain the geometry of the 
outflowing gas created during the tidal disruption event. 

In addition to emission lines from the back edge of 
the outflowing gas at late times, the photosphere of the 
super-Eddington outflow may at early times show strong 
absorption lines, particularly in the ultraviolet (much like 
the photosphere of a star); these lines would likely be highly 
blueshifted relative to the lines of the host galaxy. Finally, 
we reiterate that we expect very little narrow forbidden line 
emission from tidal disruption events: the outflowing stellar 
debris is too dense to produce forbidden lines, and there is 
insufficient time to photoionize ambient lower density gas 
far from the BH. It is possible that there is ambient gas in 
the galactic nucleus sufficiently close to the BH to produce 



forbidden lines on a <J I yr timescale, but the prevalence of 
such gas is currently poorly understood. 



6.3 Observational Considerations 

Candidate detections of tidal disruption flares have thus far 
been selected by their UV and soft X-ray emission (pre- 
dominantly via GALEX, ROSAT and XMM- Newton) . The 
emission at these wavelengths is primarily produced by 
the accretion disk, which is brighter for larger Mbh. As 
a result, these surveys are most sensitive to BHs having 
Mbh ~ fO 7 - 10 8 M Q . UV and X-ray selected events are 
likely to be discovered somewhat after the initial period of 
super-Eddington fallback, because the outflow during that 
phase probably precludes direct observation of the underly- 
ing accretion disk from many viewing angles. This outflow- 
ing material could also be a significant source of obscuration 
even at late times when the fallback rate is sub-Eddington. 

The accretion disk emission we predict for UV-selected 
events is broadly consistent with the GALEX candidates: 
UV luminosities of ~ few x 1 43 ergs" 1 , optical luminosi- 
ties of ~ few x 1 41 ergs -1 , blackbody temperatures of 
T ~ few x 1 4 — few x 10 5 K, and bolometric luminosities 
of Lboi ~ I0 45 ergs -1 . The events selected from soft X-rays 
have less data, but typically have soft X-ray luminosities 
of ~ 10 43 — I0 44 ergs -1 . This emission may be from the 
accretion disk, at energies just above the blackbody peak, 
or may be from an X-ray power-law tail with I — 10% of 
the bolometric luminosity. In addition, as discussed in H5.2I 
our model is consistent with ROSAT and GALEX rate es- 
timates, provided that the tidal disruption rate per galaxy 
for BHs with Mbh ~ IO r M - I0 8 A'/ Q is 7 ~ I0 -5 yr -1 . 
This constraint already suggests that galactic nuclei in the 
nearby universe are relatively spherical, rather than triaxial, 
because the expec ted disruption rate is significantly higher 
in the latter case (jMerritt fc Poonll2004l ). 

As described in S|4] the spectral signature of the equato- 
rial stellar debris is a transient spectrum of broad emission 
lines shifted in wavelength relative to the host galaxy. We do 
not expect forbidden lines (e.g. [Nil], [SII], [OI], [OIII]) to 
be present, because the density in the unbound material is 
too highQ For several reasons, it is not surprising that this 
spectral signature has yet to be seen. First, the tidal flare 
candidates are likely from relatively high-mass BHs; in those 
events, the unbound stellar debris subtends a small solid an- 
gle (eq. |25| i and so the emission lines should be <J 1% of 
the bulge luminosity. Future optically-selected tidal flares 
are more likely to show detectable lines. It is also impor- 
tant to note that most standard searches for AGN in opti- 
cal/infrared surveys use forbidden lines to identify nuclear 
activity, and have not specifically looked for faint, broad 
lines offset from the host galaxy's lines. Tidal disruption 
events may yet be hiding in archival spectroscopic data. 

We predict that outflows during the super-Eddington 
fallback phase have peak optical luminosities of ~ I0 43 — 
I0 44 ergs -1 and characteristic decay timescales of ~ 10 days 



7 The 5007 A line of [OIII] does appear at late times for the lowest 
Mbh an d Rp, where the densities fall below this ion's critical 
density after about a year. 
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(Figs.l4lfc lll[) . These events are sufficiently bright that a nat- 
ural concern is whether our predictions can already be ruled 
out by optical supernova searches such as the Supernova 
Legacy Survey and Stripe 82 in the Sloan Digital Sky Sur- 
vey. Although a careful search of archival data is clearly war- 
ranted, we do not believe that current observations are nec- 
essarily that constraining, for two reasons. First, the outflow 
phase can be relatively brief and many survey cadences may 
be insufficient to find these events. Most importantly, how- 
ever, tidal flares could be readily mistaken for AGN and thus 
discarded. Indeed, most supernova searches discard galactic 
nuclei in order to avoid confusion with AGN and optimize 
their probability of detecting supernovae. 

Optically detecting a tidal flare may require disen- 
tangling the flare emission from that of the BH's host 
galaxy. For example, at a distance of 300 Mpc, a ground- 
based optical survey with a resolution ~ 1" should just be 
able to resolve a kiloparsec-sized bulge. Bulges are found 
to be ~ 700 times more massive than their central BHs 
l|Haring fc Rbd [20041 ) . Super-Eddington outflows typically 
shine at 10 43 — 10 44 ergs _1 (see Fig- [4j , so this phase should 
be at least as bright as the host galaxy and readily detectable 
given sufficient attention to sources in galactic nuclei and 
careful screening to rule out an unsteady AGN. By con- 
trast, the typical optical luminosity of the accretion disk 
itself is 10 40 - few x 10 41 ergs~ 1 (see Fig. |3j. The accre- 
tion disk would thus brighten the host bulge by only a few 
percent. Photometric detections of late-time flares will re- 
quire very careful bulge subtraction. As a result, shallow, 
wide-area surveys such as PTF are more likely to find the 
late-time disk emission than narrow deep surveys such as 
the Pan-STARRS MDS. As an additional complication to 
finding tidal flares, optical extinction in galactic nuclei can 
be significant (although less than at UV or soft X-rays); as 
a result, some fraction of optical tidal flares may not be de- 
tectable due to obscuration. Infrared surveys such as SASIR, 
which are also very sensitive to tidal disruption events (Table 
[TJ, will be particularly immune to the effects of obscuration. 

Type II supernovae in the nuclear regions of galaxies 
may be confused with tidal disruption events, as both have 
quite blue colors. For sources at ~ 300 Mpc, we estimate 
that such supernovae will occur within ~ 1" of the galactic 
nucleus at a rate of ~ 10 -4 yr _1 , perhaps an order of magni- 
tude more often than tidal disruption events; at higher red- 
shift the contamination from supernovae will be more signifi- 
cant, but follow-up imaging at high spatial resolution and/or 
spectroscopic follow-up should help classify these events and 
distinguish tidal flares from nuclear supernovae. 

6.4 Astrophysical Applications 

Theoretical calculations of the tidal disruption rate per 
galaxy, 7, vary substantially and can have complicated 
dependences on BH mass and pericenter distance (e.g., 
iMagorrian fc Tremainefl999l '). The rate model we implement 
is consistent with theoretical estimates and is sufficiently 
simple that the reader can easily scale our results to different 
model parameters (eq. |32)i. We have assumed that the rate 
at which stars enter the disruption zone (_Rs < Rp < Rt) 
is independent of BH mass, and constant with ln_R p . A star 
may venture deep into the disruption zone (Rp ~ Rs) on 
its last orbit if its change in angular momentum over one 



dynamical time is large enough — at least of order the max- 
imum angular momentum for disruption. This condition is 
satisfied in the full loss cone regime, and marginally satis- 
fied in the outskirts of the diffusive regime. For realistic stel- 
lar density profiles, the disruption rate is dominated by th e 
boundary between these two regimes fe.g.. lAlexan"derll2005l ). 
so many stars probably do take large enough angular mo- 
mentum steps to arrive at R p <C Rt. In the diffusive regime, 
the disruption rate per lni? p varies weakly with \nR p , con- 
sistent with our assumption in [JS] Given the sensitivity of 
the optical- infrared emission from tidal flares to R p , upcom- 
ing surveys should significantly improve our knowledge of 
the stellar dynamics in galactic nuclei. 

Our results demonstrate that optical transient surveys 
will be quite sensitive to the lowest mass BHs in galactic 
nuclei, both because of the outflows produced when the fall- 
back rate is super-Eddington and because of the large angle 
subtended by the equatorial stellar debris. Such BHs are 
otherwise difficult to detect because their host galaxies are 
faint, it is difficult to resolve their small spheres of influence, 
and even when they are active, their Eddington luminosities 
are low. The space density of 10 5 — 1O 6 M0 BHs and the 
stellar density profiles in the galaxies they in habit are only 
mode rately well-constrained at present (e.g., I Greene fc Hoi 
l2007h . as is the role of tidal disruption in growing these BHs. 
Optical searches for tidal flares should thus prove to be a 
powerful probe of low-mass BHs and their host galaxies. 
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